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A new partial gauge fixing condition for the abelian projection is introduced. It is based on the lowest-lying 
eigenvector of a covariant Laplacian operator. This gauge is smooth and free of lattice Gribov copies. These 
properties are important for an unambiguous computation of the abelian projected gauge field configuration. 



1. INTRODUCTION 

The abelian projection [ |lj remains the most 
popular framework for lattice Monte Carlo stud- 
ies of monopoles and confinement. There are two 
uncomfortable aspects to the partial gauge fixing 
involved, though. One is the apparent preference 
for one particular gauge, the Maximally Abelian 
Gauge (MAG) [ |J, the other is the presence of 
lattice Gribov copies in this gauge. 

The situation with respect to the first point 
can be loosely summarized as follows. The MAG 
"works well", the temporal ("Polyakov") gauge 
"works sometimes", and the other gauges "do 
not work".Q (See the reviews [ @, |L p| f° r re- 
sults and interpretations.) This may have its 
origin in the different smoothness properties of 
these gauges. The MAG is a smooth gauge, the 
Polyakov gauge is smooth in the time direction 
only, while the other gauges may not be smooth 
enough. Smoothness means that the link matrices 
are relatively close to unity, which is important 
for interpreting them in terms of continuum gauge 
fields, and hence for extracting their abelian part. 

Turning attention to the apparently preferred 
MAG then, one is confronted with the second 
issue: In practice, implementation of the MAG 
is complicated by the presence of lattice Gri- 
bov copies, corresponding to different local min- 
ima of the gauge fixing functional. Unambigu- 
ous computation of the gauge fixed configuration 
is impossible, and certain quantities such as the 



*To appear in the proceedings of "Lattice 96", St. Louis, 
Missouri, USA, June 4-8, 1996. 

Work supported by EC contract ERBCHBICT941067 and 
by DGICYT project AEN 94-218. 

"However, for some new gauge conditions, see ref. [gj. 



monopole density are fairly sensitive to this am- 
biguity [ 0, H . One would like to have a smooth 
gauge without lattice Gribov problem. 

An analogous problem arises in conventional 
Landau gauge fixing. Its lattice implementation 
is smooth but there are lattice Gribov copies. To 
circumvent this problem, Vink and Wiese intro- 
duced "Laplacian gauge fixing" [ ^| , which shares 
the smoothness properties of the Landau gauge 
but avoids lattice Gribov copies. In this pro- 
posal, the gauge transformation matrices are de- 
termined in terms of the lowest-lying eigenvectors 
of the covariant Laplacian in the background of 
the given gauge field configuration. Recently this 
method was studied in practice [ lOL and its per- 
turbative formulation was given [ . 

Here I take this idea over to the abelian projec- 
tion, constructing a procedure for partial gauge 
fixing which is smooth and free from ambigui- 
ties. Section 2 introduces the gauge, discusses its 
expected merits and briefly describes its pertur- 
bative continuum formulation. Numerical results 
are presented in section 3. 

2. LAPLACIAN ABELIAN GAUGE 

To introduce the Laplacian Abelian Gauge 
(LAG) it is convenient to start from the MAG 
in its spin- model formulation [ I will limit 
myself to the SU(2) case here. 

The Maximally Abelian Gauge [ ||] for a link 
configuration {U^ yX \ is defined as the configura- 
tion {Ufj,, x = OxU^^fl^^} where minimizes 
the functional 



Su(Q) = 



E{l-|Tr 



}, (1) 
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which can be written in the form 

M<t>) = E i 1 - ^ Tr ■„():.,:} (2) 

= E^-EW^+A}. (3) 



with the definitions 



(4) 
(5) 



is the link matrix in the adjoint representa- 
tion and 4>x S SU(2)/U{1) ~ 5 2 is a three-vector 
of unit length, neutral under abelian gauge trans- 
formations, as expected. The gauge fixing func- 
tional (^) is nothing but the latticized covariant 
kinetic action J i(D M 0) 2 of a spin field <f> in the 
background of the given gauge field. The spin 
configuration which minimizes this action deter- 
mines the MAG. Finding this absolute minimum 
is difficult because of the length-one constraints 
on the individual spins. The usual iterative local 
algorithms often get stuck in a local minimum. 
This means that the result of the gauge fixing de- 
pends on the particular algorithm used and on 
the starting point on the gauge orbit: the gauge 
fixing procedure is ambiguous. 

The idea of Laplacian gauge fixing is to: 

1 . Minimize Su (4>) without taking into account 
the constraints \\4> x \\ = 1- This amounts to find- 
ing the eigenvector belonging to the lowest 
eigenvalue A of the covariant Laplacian — D x t(R)' 

2. Write the solution as <^>" = p x <j)% and take 
for the gauge transformation: <& x —^'t'x^a = 

Cl+asClx (of course Cl x is determined up to the 
residual U(l) freedom only). 

This procedure is unambiguous because the 
computation of eigenvectors can be done to the 
precision required, and because the procedure is 
gauge covariant by construction: under a gauge 
transformation V of the starting configuration 
U, the Laplacian operator — 0(R(U)) and its 
eigenvectors transform accordingly, such that the 
gauge fixed configuration U is unchanged (again, 
up to residual abelian gauge transformations). 

The procedure is ill-defined only if the lowest 
two eigenvalues of the covariant Laplacian coin- 
cide, or if <h x = for some x. However, the set 



of such configurations has measure zero. In prac- 
tice, we might get into problems when either the 
difference between the two lowest eigenvalues, or 
4> x for some x, is zero within the numerical pre- 
cision of the computer, but, as it turns out, this 
never occurs. 

The possibility of a zero in <fi (not necessar- 
ily at a lattice site) is actually quite interesting. 
It means that the gauge fixing is ill-defined at 
that point, which is precisely what identifies a 
magnetic monopole in the abelian projection. In 
fact, the 't Hooft-Polyakov monopole (dyon) in 
the radial gauge satisfies the continuum equiva- 
lent (see below) of the LAG, with the solution for 
4> equal to the Higgs field (A4) of this configura- 
tion (which has a zero at the origin)! 

In the continuum limit, the stationarity condi- 
tions corresponding to the LAG are 



E 



-dl + {Alf + {Alf)p = Xp, 



E ^ T iA l 



(G) 



(7) 



(Note that p(x) depends on the gauge orbit as a 
whole only.) This corresponds to minimization of 
the quantity 



J V P 2 [(Alf + {Alf] 
IvP 2 



(8) 



(its minimum value being the smallest eigenvalue 
A of the covariant Laplacian) or, on the lattice, 



1 - P ({U}) p ({u)) i Tr 



}-(9) 



Note that in all these formulas one recovers 
the corresponding expression for the MAG [ [ll| 
by setting p equal to 1. In fact, in the contin- 
uum limit (3 — » 00, the LAG- fixed configuration 
is characterized by U^x — > 1 and /5, -» 1, and 
one sees that LAG and MAG converge to each 
other. 

At finite /?, the function p acts like a kind of 
"measure of local smoothness" : eq. (||) suggests 
that p(x) will be small when (A l ^) 2 (x) + (A 2 ) 2 (x) 
wants to be large. For example, in the centre of 
the 't Hooft-Polyakov dyon the non-abelian field 
components A 1 ^ 2 blow up while p — 0. Such di- 
vergences thus bear a relatively low penalty in 



the gauge fixing functional (||): the requirement 
of smoothness appears to tell the LAG to treat 
monopoles more mildly than the MAG does. 



3. NUMERICAL TESTS 

The following table contains some results for a 
6 4 lattice at various (3 values. 



p 


2.2 


2.4 


2.6 


2.8 


Ai 


2.123(3) 


1.872(8) 


1.677(2) 


1.529(2) 


A 2 


2.154(3) 


1.926(8) 


1.759(2) 


1.631(3) 


Ap 


0.512(3) 


0.402(1) 


0.293(2) 


0.243(2) 


#L 


492(6) 


149(13) 


28.2(8) 


6.6(4) 


#M 


428(4) 


103(10) 


16.8(5) 


4.2(3) 



Ai ; 2 are the two lowest eigenvalues of the co- 
variant Laplacian, and Ap is the average fluctu- 
ation of p over a configuration. For (3 — > oo, Ai 
tends to zero, A2 to one (becoming 24-fold de- 
generate), and Ap to zero. However, one should 
keep in mind that this limit is unphysical since it 
corresponds to zero physical volume. 

#l and #m ar e the numbers of dual links 
transmitting monopole current, for LAG and 
MAG respectively. Note that the monopole num- 
bers are somewhat higher in the LAG. This 
should not be interpreted as a sign that the LAG 
is "worse" than the MAG, as might be suggested 
by the empirical fact that "bad" gauges usually 
have very high monopole numbers associated to 
them. Rather, it may be regarded as an indi- 
cation of the number of monopoles present in a 
situation of optimal smoothness. 

Fig. 1 shows a comparison of the monopole lo- 
cations as determined by the LAG and the MAG, 
for a randomly chosen configuration at (3 — 2.4. 
It is interesting to see the large overlap, which is 
a sign of the similarity between the two gauges. 
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Figure 1. Monopole loops for a 6 4 configuration 
at P = 2.4, using LAG (a) and MAG (b). Open 
circles mark sites in the boundary of the (dual) 
lattice, full circles indicate sites in the interior. 
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